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Abstract 



We report a simple method to accurately determine the thresh- 
old and the exponent v of the Bak-Sneppen model and also inves- 
tigate the BS universality class. For the random-neighbor version 
of the BS model, we find the threshold x* = 0.33332(3), in agree- 
ment with the exact result x* = 1/3 given by mean- field theory. 
For the one-dimensional original model, we find x* = 0.6672(2) 
in good agreement with the results reported in the literature; for 
the anisotropic BS model we obtain x* = 0.7240(1). We study the 
finite size effect x*{L) — x*{L oo) oc L^"^, observed in a sys- 
tem with L sites, and find = 1.00(1) for the random-neighbor 
version, z/ = 1.40(1) for the original model, and = 1.58(1) for 
the anisotropic case. Finally, we discuss the effect of defining the 
extremal site as the one which minimizes a general function /(x), 
instead of simply f{x) = x as in the original updating rule. We 
emphasize that models with extremal dynamics have singular sta- 
tionary probability distributions p{x). Our simulations indicate 
the existence of two symmetry-based universality classes. 
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I. INTRODUCTION 



Power-law distributions are often found in natural and social systems. Ex- 
amples include the distributions of rain event sizes and durations [1-4] , earth- 
quake sizes [5,6], internet connections [7], and food webs [8]. Why do such 
diverse phenomena exhibit scale invariancc, as reflected in power-law distri- 
butions? A concept introduced by Bak and collaborators to account in part 
for scale-invariance in nature is self-organized criticality (SOC), which asserts 
that "systems that are far from equilibrium evolve to a dynamical attractor 
poised at criticahty" [9]. 

Over the last two decades, many models exhibing self- organized criticality 
have been proposed: sandpile models [6,10,11], the OFC model for earthquakes 
[12], forest fire models [13], the ballistic deposition model of rough surfaces [14], 
and the Bak-Sneppen evolution model [15]. Many SOC models are defined in 
such a manner that they are trapped at criticality [16]. In particular, sandpile 
models are maintained at the critical point by "slow driving in a system ex- 
hibiting an absorbing-state phase transition with a conserved density" [16]. A 
different path to SOC is so-called extremal dynamics. In this case, the system 
is critical provided that an 'extremal condition' is respected [16-18]. 

The Bak-Sneppen evolution model is one of the simplest models exhibiting 
SOC via extremal dynamics. For this reason, it has atracted much attention 
in the statistical physics community and has been studied through various 
approaches, including simulation [19,21,22], theoretical analysis [23-25], prob- 
abilistic analysis (run time statistics) [26,27], renormalization group [28,29], 
field theory [30] and mean- field theory [16-18,31]. Some variants have been 
proposed, for example the anisotropic BS model [32,33]. 

In this paper, we investigate the scaling behavior of the Bak-Sneppen 
model. We study the finite size scaling = kL , where xj^ is the 

threshold for a one-dimensional system of L sites and /c is a constant. This 
provides an accurate estimate of the threshold and of the exponent u. We 
also discuss which modifications of the updating rule preserve the universality 
class. In agreement with Ref. [18], we find that extremal dynamics leads to 
a singular stationary distribution p{x), and that models with symmetric up- 
dating rules belong to the same universality class of the original model. The 
specific form of the distribution p{x) is, however, nonuniversal. 

The balance of this paper is organized as follows. In Sec. II, we define the 
BS model along with the random-neighbor and the anisotropic BS versions. We 
define also a "cosine BS model" . We calculate the threshold and the exponent 
u of the first three models in Sec. III. In Sec. IV, we study the BS universality 
class. We summarize our conclusions in Sec. V. 

II. MODELS 

The Bak-Sneppen model [15] is defined on a d-dimensional lattice with L'^ 
sites and periodic boundaries. Each site bears a real- valued variable Xj e [0, 1). 
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At time zero, all the Xi are assigned independent random values drawn from 
a distribution uniform on [0, 1). At each step of the evolution, the site x^, 
that is the global minimum of the {xi}, along with its 2d nearest neighbors, 
are replaced by new random numbers taken from the same distribution. This 
process is repeated indefinitely. 

In the random- neighbor version of the model, all sites are considered neigh- 
bors and, at each time step, the site m along with two randomly chosen neigh- 
bors are substituted with new random numbers. All the other definitions 
remain the same. A further variant of the model is the anisotropic BS model. 
In this case, only the site m and its right neighbor are replaced at each time 
step. We also study a modified Bak-Sneppen model, that differs from the 
original in the way the site Xm is identified. We baptize this the "cosine BS 
model" . Instead of minimizing f{x) = x as in the original model, x^ minimizes 
f{x) = cos(2n7ra;) for some integer n, so that / has n minima on the interval 
[0,1]. Finally, in order to investigate the effect of symmetry in the updating 
rule, we define an "anisotropic cosine BS model" , in which we have the spatial 
anisotropy of the anisotropic BS model and the extremal site defined as in the 
cosine BS model. 



Simulation of the BS model shows that, after a transient time, the system 
achieves a statistically stationary state in which the density of Xi vanishes for 
X less than a certain threshold x* and is uniform above x*. In other words, 
for an infinite system and in the stationary state, the distribution p{x) is the 
step-fuction 



This singular probability density is associated with scaling properties rem- 
iniscent of a critical point; the upper critical dimension rfc = 4 [22]. The 
Bak-Sneppen model does not appear to belong to any of the known universal- 
ity classes of critical phenomena, and may be said to define its proper class. 

Accurate determinations of x* were reported by Grassberger [19] and 
Paczuski and coworkers [9], who found x* = 0.66702(8) and x* = 0.66702(3), 
respectively. We introduce here a third method to determine x*. For a finite 
system (see Fig. 1), the distribution p{x) is continuous, i.e., the step func- 
tion singularity is rounded [20]. Nevertheless, as is apparent from the figure, 
p{x) = Cl, constant, for x well above the transition region. If p{x) is a step 
function, normahzation imphes that C{1 — x*) — 1. This suggests that we 
define a finite-system threshold via x}^ = 1 — 1/Cl- If we now assume that the 
Bak-Sneppen model obeys the finite-size scaling relation 



III. DETERMINING THE THRESHOLD 




X < X* 
X > X* 



(1) 



* * 7 7- — l/l/ 



(2) 
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where A; is a constant and x^^ is the threshold for an infinite system, we have 
a simple method for determining the threshold of Bak-Sneppen-like models. 

We simulate the Bak-Sneppen model, the anisotropic BS model, and the 
random-neighbor version on lattices of L = 63, 125, 250, 500, 1000, 2000, 4000 
and 8000 sites (see Figure 1). We estimate the probabihty density p{x) on the 
basis of a histogram of barrier frequencies, dividing [0,1] into 100 subintcrvals. 
Histograms are accumulated after Ngt time steps, as required for the system to 
relax to the stationary state. (We find that for the original and the anisotropic 
models Ngt — 10^ sufficient to reach the stationary states for even the largest 
system studied, while for the random-neighbor version Nst — 10^ is sufficient.) 
We make the histogram every L time steps in runs of 10^ steps total. 

To test the accuracy of the method described above, we applied it to the 
random- neighbor version of the model, for which the exact result x* = 1/3 
is known from mean-field theory [16-18,31]. Analyzing our data for x*j^ vs. 
L using Eq. (2), we find x* = 0.33332(3) and ly = 1.00(1). This results 
agrees with the mean-field theory prediction and validates our method. For the 
original BS model, we find x* = 0.6672(2) and u = 1.40(1). (The power law fit 
is shown in Fig. 2.) In accord with this are the results x* = 0.66702(8) reported 
by Grassberger [19], x* = 0.66702(3) found by Paczuski and collaborators [9], 
and u — 1.36 [34]. Finally, for the anisotropic BS model, we obtain x* — 
0.7240(1) and v — 1.58(1). To our knowledge, this is the first result for the 
threshold and the exponent i/ of the anisotropic BS model. 

IV. THE BS UNIVERSALITY CLASS 

As noted in the Introduction, extremal dynamics is one path to a self- 
organized critical state. Two questions therefore arise: 1) What are the trade- 
marks of extremal dynamics? and 2) What other models belong to the uni- 
versality class of the Bak-Sneppen model? In order to answer these questions, 
Garcia and Dickman [18] studied the effects of updating via x[ = xf (with 
a = 1/2 or 2), instead of replacing these variables with random numbers. They 
found that these updating rules lead to quite distinct stationary probability 
densities p{x), which nonetheless exhibit one or more discontinuities. Another 
caracteristic shared by all the variants studied in [18] is that the extremal site 
always belongs to the "prohibited region" for which p{x) = 0. Although p{x) 
and consequently the threshold are nonuniversal, all variants studied belong 
to the universality class of the original model. Moreover, it has been shown 
[33,35,36] that models that break the symmetry of the original model fall into 
the universality class of the anisotropic BS model. 

A related question is the robustness of BS scaling under modifications of 
extremal dynamics. It is known that scale invariance is lost when the ex- 
tremal condition is relaxed, for example by introducing a nonzero probability 
for choosing a nonextremal site to be updated [16-18]. In this context it is 
interesting to note that updating the M smallest barriers (with M > 1), at 
each step changes the scaling properties drastically [37]. 
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As representatives of a further broad class of variants of the Bak-Sneppen 
model, we study here the isotropic and the anisotropic cosine BS models pro- 
posed in Sec. 11. We simulate these models (with n = 5) on a lattice of 
L — 2000 sites and find the distribution p{x) shown in figure 3 for the isotropic 
case. (In the anisotropic version, a similar distribution is found.) It is natural 
to regard the rounding of the singularities as a finite-size effect. Thus we find 
evidence that these updating rules also lead to singular probabilty densities. 
We have also investigated the distribution Pj{r) that sucessive updated sites 
are separated by a distance r. This obeys the power law Pj{r) ~ r~'^ with 
TT = 3.23(2) in the original model [32], and tt = 2.401(2) in the anisotropic BS 
model [33]. We performed simulations of the cosine BS models using 10^ time 
steps, yielding tt = 3.2(1) for the symmetric updating rule and vr = 2.38(3) for 
asymmetric updating rule. Therefore a dynamics that respects the symmetry 
of the original model (that is, spatial isotropy), and that does not introduce 
any new conserved quantity, leads to the same exponents as the original model. 
On the other hand, asymmetric dynamics leads to the anisotropic BS model 
universality class. 



V. CONCLUSIONS 

We verify the finite-size scaling relation a;^ — x*^ = kL'" in the Bak- 
Sneppen model, and use it to obtain accurate values for the threshold and for 
the exponent of the original, anisotropic and random-neighbor versions. We 
also study the effect of defining the extremal site as the one that minimizes 
a function f{x) with multiple minima. Simulations show that this does not 
affect the universality class and strengthens the result of Ref. [18] that extremal 
dynamics leads to singular probability densities. Finally, our results indicate 
that variants that preserve the spatial isotropy of the original model belong to 
the same universality class, while models which break that symmetry fall into 
the anisotropic BS model universality class. 
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FIG. 1. Stationary density p{x) for the original Bak-Sneppen model in one di- 
mension, for various system sizes. 
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FIG. 2. Power law fit for (x^ — vs. 1/L for the original BS model. The 
slope of the straight line is 1.40. 







P(x) 








g(x) 


=f(x)+4=cos(10Hx)+4 

\ /' '\ \ /' 




V/' 




\ 


ill! ill 

1 , 1 , 1 , 1 , 1 , 1 , i 



FIG. 3. Solid line: stationary density p{x) for the cosine BS model (n = 5) for a 
system of L = 2000 sites. Dashed line: function f{x) = cos(107rx) that is minimized 
in the cosine BS model, (/(x) is shifted for better visualization.) For further details, 
see Section II. 
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